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We address ultracold A'^-atom degenerate cascade systems and show that stationary subradiant 
states, already observed in the semiclassical regime, also exist in a fully quantum regime and for 
a small number of atoms. We explicitly evaluate the amount of stationary entanglement for the 
two-atom configuration and show full inseparability for the three-atom case. We also show that 
a continuous variable description of the systems is not suitable to detect entanglement due to the 
nonGaussianity of subradiant states. 
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I. INTRODUCTION 

Atomic subradiant states have recently gained wide at- 
tention because of their exceptionally slow decoherence 
[l|,|2|. This stability of quantum superpositions inside the 
subradiant subspaces originates from the low probability 
of photon emission, which means very weak interaction 
between the atoms and their environment. Hence, the 
subradiant states span a decoherence-free subspace [3|-|5| 
of the atomic Hilbcrt space, and consequently, can be- 
come important from the viewpoint of quantum compu- 
tation [aQ. 

Subradiance in an extended pencil-shaped sample has 
been observed in an unique experiment in 1985 by 
Pavolini et al. [SJ. Among different schemes of multi- 
level systems, Crubellier et al.[9 12] investigated in de- 
tail the arising of subradiance in three-level atoms with 
a degenerate transition. More precisely, they considered 
the case of three-level atoms with two transitions shar- 
ing either the upper level ('A configuration) or the lower 
level ('F' configuration) or the intermediate one ('cas- 
cade' configuration), when both the transitions have the 
same frequency and the same polarization. The first two 
configurations may be experimentally realized by rather 
specific atomic level systems, having a small hyperfine 
structure, whereas the degenerate cascade configuration, 
which contains two cascading transitions of the same fre- 
quency, would in fact be encountered in atomic-level sys- 
tems only in the presence of an external field which suit- 
ably modifies the atomic frequencies jlfll- Recently, the 
same degenerate cascade among three equally spaced lev- 
els has been investigated for collective momentum states 
of an ultracold atomic gas in a high-finesse ring cavity 
driven by a two- frequency laser field 13*. This system is 
particularly attractive for the observation of subradiance, 
since the absence of Doppler broadening and collisions at 
sub-recoil temperatures avoid other undesirable decoher- 
ence mechanisms. Furthermore, the experimental control 
of the external parameters, as for instance the intensities 
of the pump fields, allows for a continuous tuning through 



different subradiant states of the system, which instead 
should result difficult or impossible for internal atomic 
transitions. 

More specifically, the system described in |13j con- 
sists in N two-level ultracold bosonic atoms (at T « 0) 
placed in a ring cavity with linewidth k and driven by 
two pump fields Eq and Ei of frequency lo and w -I- A, 
respectively. The atoms scatter the photons of the two- 
frequency pump into a single counterpropagating cavity 
mode. The pump frequency w is sufficiently detuned from 
the atomic resonance wq to neglect absorption. The scat- 
tering process can be described by two steps. In the first 
step the atom, initially at rest, scatters the lo pump pho- 
ton into a photon with frequency Usi = cj — w,., and 
recoils with a momentum p = hq (where q — k — kg 
is the transferred momentum, lOj. = q^ jIM is the re- 
coil frequency and M is the atomic mass). In the sec- 
ond step the atom scatters the w -I- A pump photon into 
a photon with frequency Ws2 = w + A — Sw^, chang- 
ing momentum from p = hq into p ' — Ihq. The recoil 
shift of SciJr arises from the kinetic energy conservation, 
i.e. b.E = (p'2 - p^)/2M = ihwr. If the frequency 
difference between the two pump fields is A = 20;^, 
then LOsi = uJs2 = ^s and a degenerate cascade between 
three momentum levels, p — m{hq), with m — 0,1,2, 
is realized. Furthermore, if the cavity linewidth k is 
smaller that 2iJr, the cavity will support only a single fre- 
quency uJs , avoiding further higher-order scattering with 
frequency lOs ± mojr- In this way the transitions are re- 
stricted to only the first three momentum states, forming 
a three-level degenerate cascade. Although the condition 
K <^ 2ujr is experimentally demanding, requiring a typi- 
cal cavity finesse larger than 10^, it is not very far from 
the state-of-art optical cavities [ij, \l^ . An important 
feature of the momentum three-level cascade is that the 
ratio of the two transition rates is e = \Ei/Eo\, i.e. it is 
an experimentally tunable parameter, contrarily to the 
case of atomic transitions where it is fixed by the branch- 
ing ratios. 

In this paper we address A^-atom degenerate cascade 



systems and show that stationary subradiant states also 
exist in a fuUy quantum regime. ActuaUy, the action of 
the coupHng to the environment drives our system into 
an entangled state without the need of isolating it from 
the rest of the universe. We solve explicitly the Mas- 
ter equation governing the system for a small number of 
atoms and evaluate the amount of stationary entangle- 
ment [16{ . In turn, our system represents an example 
of robust generation of steady state entanglement among 
atoms via the interaction with the electromagnetic field 
which act as an engineered environment [ITf . The paper 
is structured as follows. In the next Section we give a 
full quantum description of the dynamics and we numer- 
ically solve a quantum master equation for N = 2(3)- 
atom systems. In the Section ITTll we study the entangle- 
ment properties of subradiant states looking at the sys- 
tem from both discrete- variable and continuous- variable 
(CV) point of view and show that CV criteria are not 
able to detect entanglement due to nonGaussianity of 
subradiant states. Section |TV] closes the paper with some 
concluding remarks. 



|sr)o = |0, 0, N) for p =: 0. If A^ > 1 these states exist iff 
0<e<l + V2. 

In order to study the time-evolution of g accounting 
for scattered photons dissipation too, in a pure quantum 
mechanical framework, we introduce the following T = 
Lindblad master equation 



dt 
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[Hint.g] + 2KC[a]Q 



(3) 



where /C[a]£i = aga* — a*ag/2 — ga*a/2 is the Lind- 
blad operator describing the cavity damping of the mode 
a. In the next subsections we solve numerically Eq.Q 
for two and three atoms in terms of total Fock states 
|no, ni, 712, n), where n is the number of scattered pho- 
tons. We will find that for every e and for k 7^ Eq.Q 
always admits a stationary solution of the form 



Qs 



Pi|sr)n(sr| + (l-Pi)|sr)oo(sr| ® |0)(0| (4) 



where |sr)o,|sr)i are the only two populated states of 
two-atom and three-atom systems. Pi = Tr {g\sr) 1 1 (sr|} 
and |0) is the vacuum photon state. 



II. DISSIPATIVE DYNAMICS AND N-ATOM 
PURE SUBRADIANT STATES 



A. N = 2 system 



The three-level degenerate cascade interaction is de- 
scribed by the following Hamiltonian operator in inter- 
action picture 



Hi, 



-ihg 



acicj 



eac2C{ — h.c. 



(1) 



for the three atomic bosonic operators, c™, with m — 
0, 1, 2 and [c„i, cjj — Sm,n, and for the single cavity radi- 
ation mode, a, with [a, a^] = 1; g is the coupling constant 
and e is the relative rate of the second transition, from 
m = 1 to TO = 2. The state of the total system is de- 
scribed by a density operator g which is defined on the 
overall Hilbert space Ti = 'Hq®'Hi®'H2® "H/, Tii being 
the Hilbert-Fock space associated to the m = i atomic 
mode and Hf the Hilbert-Fock space associated to the 
radiation mode a. 

Introducing the Fock states |no,rii,ri2), where rii is the 
number of quanta in the m — i mode, it has been demon- 
strated in [13, [13 that such a A^-atom system, with N 
even, admits A'^/2 subradiant pure states. 
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where 2Fi{a,b;c; z) is the hypergeometric function and 
p = 1, . . . , A^/2, at whom we can add the ground state 



From Eq. ^ it is straightforward to check that the two- 
atom system admits the subradiant state 

|sr)i = -=i==(|0,2,0)-V2e|l,0,l)). (5) 

Vl + 2e^ 

We choose as initial condition for Eq.dS]) gi — 
|2,0,0,0)(2,0,0,0|. FigUKtop) shows the time evolution 
of the populations (Ni) = Tr[c]ci£i] (with i = 0, 1, 2) and 
(N) — Ti[a''ag] for chosen values of k and e: we observe 
that (A'o) and (A^i) tend to non-zero values after the pho- 
tons have escaped from the cavity. This result suggests 
that the atoms do not decay to the lower momentum state 
TO, = 2, but they arrange themselves in a superposition 
of the three momentum states. Figl5] (top) shows the 
time evolution of the probabilities Pi — TT[g\sr)ii{sr\], 
for j = 0, 1, for the same values of k and e. Asymptoti- 
cally, Pq + Pi = 1 and g approaches the steady-state den- 
sity matrix gg of Eq. ^ , where Pi in general is a function 
of the parameters e and k. A similar behavior may be 
observed for any choice of the interaction parameters. In 
the bad-cavity limit, k,^ g, Eq. ([3]) is well approximated 
by the following effective superradiant master equation 



—— = Vb gaS^ 
dt 



-{S+S^ga 



QaS^S' 



(6) 



where ga is the atomic density operator, S — c[co + 
ec^ci and F = g^ / k. Solving this equation analytically 
we get an expression for the time-evolution of Pi and in 
the stationary limit we find a F-independent expression 

2{l-e^f 



9e2 -H 2(1 - e2)2 



(7) 



This expression is in excellent agreement with the solu- 
tions of Eq. dSl) for any possible value of e in the bad 
cavity regime. In Fig|3] (top) we see that Pi decreases 
from unit as e increases and it vanishes for e = 1, for 
which the degenerate three-level cascade has equal tran- 
sition rates and the atoms decay superradiantly to the 
ground state. For e > 1 the subradiant probability in- 
creases and reaches asymptotically the unity. For large 
values of e the steady-state is no longer a mixed state, 
e,«|l,0,l)(l,0,l|. 




decreases increasing e and vanishes for e = 1, for A'' = 3 
Pi increases, it reaches a maximum around e « 0.5 and 
decreases again falling off to zero for e = 1 . For e > 1 the 
behavior of the two plots is quite similar. 
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FIG. 2. Time evolution of Pq (blue dotted) and Pi (red) with 
t in unit of 1/g. Left: N ^ 2, k = 0.2g, t = 0.3; Right: iV = 3, 
K = 0.3(7, £ ~ 0.5. 
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FIG. 1. Time evolution of (iV,) for m = (red), m = 1 (blue 
dashed), m — 2 (green dotted) and {Nph) (black dot-dashed) 
with t in unit of 1/g. Left: N = 2, k = 0.2g, e = 0.3. Right: 
Af = 3, K = 0.3g, 6 = 0.5. 



C. Large TV and p-degeneracy 

Let us consider now the A^ ^ 1. In [1J| the following 
relation linking the p-index and e has been given: 



B. 



A'^ = 3 system 



In this case we cannot directly use the eq.([2|) which 
applies only for N even. The equation S~\^j) = 0, where 
\'ip) is a generic three-atom state, is satisfied by the sub- 
radiant state 

■ ' ^ =(|0,2,l)-2e|l,0,2)) (8) 



\sr I 



Vl + 4e2 



FigH] (bottom) shows the time evolution of the popu- 
lations {Ni) with i = 0,1,2 and (A^) for fixed values 
of e and k, obtained solving Eq.Q with initial condi- 
tion ft = |3,0,0,0)(3,0,0,0|. As in the A^ = 2 case the 
mean occupation number of every atomic mode saturates 
to non-zero values and, asymptotically, the only atomic 
states surviving are the subradiant state \sr)i and the 
ground state |sr)o = 1 0,0,3). Hence also in this case we 
can write the reduced steady-state density operator in the 
form (g]). Fig 12] (bottom) shows Pi vs t and Fig (3] (top) vs 
e for a particular choice of k. There are some important 
differences and analogies with the N — 2 case. First we 
note that this time Pi vanishes for e = 0, whereas in the 
N = 2 case it has a maximum and is equal to unity. This 
is a consequence of the fact that the two-atom subradiant 
state ([5]) contains the state 1 0,2,0) which is the ground 
state for e = (only the m — -^ m ~ 1 transition can 
take place). On the contrary, for A^ = 3 the steady-state 
density operator does not contain the state |0, 3, 0) which 
is the three-atom analogue; moreover, it is a linear super- 
position of states with n2 ^ so, if the second transition 
is forbidden, it is impossible to populate these states ob- 
serving subradiance. Secondly, whereas for N = 2 Pi 



N /l-2e" 



^ = Y 



p 



N 



,2N 2 



1-e' 



< e < I/VS (9) 

l/%/3<e<l + V2 (10) 



Let us pay more attention to formulas pO)) . For l/VS < 
e < \/3 it is easy to prove that they exist two different 
values of e corresponding to the same value of p (except 
for e = 1 yielding p = i.e. the ground state). This 
means that, once p is fixed, the two subradiant states 
l'5'')p(co)' k^)p(ei) ^^6 both acceptable, provided that they 
are actually different states. Now, let us assume that, for 
a fixed p, it is possible to build the following two states: 



1 



\Wp ') = -^i\sr)p{eo) ± \sr)p(e,))- 



(11) 



It is straightforward to show that the two following states 
are orthonormal 



14^^) = 



1 



VliS 



1^;,^ 



(±) 



(12) 



We 



where a = p{eo){sr\sr)p(,^) = pie,){sr\sr)p(,^) G 

can compute the total kinetic energy of the condensate 

and also the energy difference between these two states 



-\j^(+)i 



?(-) 



Apl = \E:y'{e) - E^p '{e)\ which turns out to be 



|A„ 



4{2N-2p-2[k{eoip)) + k{ei{p))]) 



V2(l + «) 
2meo{p)) + Heiip))] 



\/2(l-a) 



(13) 



where 



fc(e,(p)) = 



p\{N -p-iy. 

2Fi(l-p,3/2;l+p-JV;e^T 
2F,i~p,l/2;p-N;ef) 



Since this parameter depends on e, in principle it should 
be possible to distinguish between the two states (IT^ . 
More interesting, by applying an external potential cou- 
pling these two states, we could realize a fictitious two- 
level system where the many-body states \(f>p ) play the 
role of the logical-basis states |0), |1). A similar system 
has been proposed in [19i | for a A configuration. 





FIG. 3. Left: steady state probability as a function of e for 
N = 2 (red) and A'' = 3 (blue dashed). Right: negative 
eigenvalue of the partially transposed state as a function of e 
for N = 2 (red) and A'^ = 3 (blue dashed). In the A'^ = 3 case 
we have set k = O.Sg. 



III. ENTANGLEMENT PROPERTIES 

In this section we address the entanglement properties 
of the state ^ cither regarded as made of three bosonic 
modes or, upon considering the selection rule induced 
by the constraint on the number of atoms, as a finite 
dimensional one. Generally speaking, there are several 
classes of entanglement in a tripartite quantum system, 
from genuine tripartite entanglement to full separability 
passing through bipartite separability, which are defined 
on the basis of all possible groupings of the three parties 
[20|. Our primary goal is to establish whether the sta- 
tionary state is fully inseparable and if it is possible to 
quantify the amount of entanglement. 

At first, we notice that our system consists of three 
atomic modes for which, being the atomic ensemble made 
of a finite and fixed number of atoms, only a finite number 
of quanta per mode is permitted from a minimum of rii — 
up to a maximum of n.i = N . This physical constraint 
corresponds to a geometrical constraint which identifies 
a subspace S of Tia — 'Hq®'Hi®'H2 defined as 



for discrete- variable quantum systems. In fact, PPT con- 
dition is sufficient to establish entanglement of a state 
and, as we will see, it is enough to show that Qs is gen- 
uinely tripartite entangled. Partially transposing with re- 
spect to a given party corresponds to a specific grouping 
of two parties with respect to the third: for example when 
we partial-transpose with respect to the m = mode we 
are grouping together the to = 1 and m = 2 modes with 
respect to ra — 0. Suppose we have partial-transposed 
with respect to jti = 0, obtaining q^° and let the eigen- 
values of g^° be {A^ '}. If it exists at least one k such 
that Afc < wc can conclude that Qs ^ X) i Pj dj ® qY ■ If 
the positivity condition is violated for q^^ and g^"^ too we 
can state that Qs is genuinely tripartite entangled. We 
start with the N — 2 case and we take the reduced den- 
sity operator Trp^ [qs\ where Qs is given by eq. ^ with Pi 
given by eq.©. The S subspace is spanned by six Pock 
states: |2, 0, 0), |0, 2, 0), |0,0, 2), |1, 1, 0), |1,0, 1), |0, 1, 1). 
If we compute the eigenvalues of every Qs' with j = 0, 1, 2 
we find that it always exists a negative eigenvalue ^2(e) 



A2{e) = - 



V2 e{e^ - 1)^ 
2 (1+ £2)2(2 + e2) 

£^(26^ + 3)^ + 2 
4V2(e2-l)2 



{No){Ni 



(14) 



where the expression on the second line provides a link 
between the negative eigenvalue, which directly quanti- 
fies the amount of entanglement, and a couple of mo- 
ments which are, at least in principle, amenable to direct 
measurements. 

In the A^ = 3 case S is spanned by the following ten 
states: |3,0,0), |0,3,0), |0,0,3), |2,1,0), |2,0,1), |0,1,2), 
|0,2, 1), |1,0,2), |1,2,0), |1,1,1) and the PPT condition 
is still quite simple to employ to establish tripartite en- 
tanglement. On the other hand, for A^ = 3 we have 
not been able to give an analytical expression for Pi and 
therefore we resort numerical evaluation of the eigenval- 
ues. As in the previous case we find a negative eigenvalue 
^3(e) which turns out to be always the same for every 
possible partial transposition and for any value of e. This 
means that also the three-atom g^ is genuinely tripartite 
entangled. Fig [3] (bottom) shows the behavior of ^jv(e) 
both for A = 2 and A^ = 3. The behavior is qualitatively 
the same for both cases. The eigenvalues decrease from 
zero and reach a minimum after which they start to in- 
crease. As we should expect they vanish for e = 1 and 
after this they start to decrease again. 



A. Continuous variable description 



S^lgen :^(A,) 



A 



i=0 



Hence, being S a finite dimension subspace, we can make 
use of the PPT {Positive Partial Transposition) criterion 



Let us now investigate the entanglement properties of 
our system as a continuous variable one [2l|, i.e. tak- 
ing into account that we are dealing with three bosonic 
modes describing the number of atoms within a given 
state (for instance momentum state). To this aim we 
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FIG. 4. Up: 5[q\ for TV = 50 vs. e accounting for relations 
®-(Il0l). Bottom: 5[q\ for iV = 2 (red) and Af = 3 (blue 
dashed); in this second case we have set k = O.Sg. 



start by reviewing few facts about the description of a 
multimode continuous variable system [23|. Upon intro- 
ducing the canonical operators qj = {cj + cA/^j^ and 

Pj = i(Cj — Cj)/^/2 in terms of the single-mode opera- 
tors Cj, j = 0, 1, 2 we may write the vector of operators 
R= (gi , pi , ^2 , P2 , <Z3 , P3 )^ ■ The so-called vector of mean 
values X ~ X[g] and covariance matrix (CM) cr = (T[g] 
of a quantum state g have components Xj — (Rj) and 
Rajk = ^{{Rj,Rk}) - {Rj){R3) where (O) = TiigO] 
and {A, B} = AB + BA denotes the anticommuta- 
tor. The characteristic function x[e](A) = Tr[gi£'(A)] 
is defined in terms of the multimode displacement op- 
erator D{X) = «'"^i£'(Aj), with A = (Ai,--- ,A„)^, 

Xj G C and D{\j) ~ expjAjcJ — A*Cj}. A quantum 
state go is referred to as a Gaussian state if its char- 
acteristic function has the Gaussian form x[qg\{-^) = 
exp<— iA crA + X J7A > where A is the real vector 

A = ($HeAi,nmAi, • • • ,9^eA„,3mA„) and Jl is the sym- 
plectic matrix Vt = i ©7=i <^2, o'2 being the y Pauli ma- 
trix. Gaussian states are thus completely characterized 
by their mean values vector and their covariance ma- 
trix. For bipartite Gaussian states, the PPT condition 
is necessary and sufficient for separability [24|, and it 
may be rewritten in terms of the CM, which uniquely 
determines the entanglement properties of the state un- 
der investigation. A mode of a Gaussian state is sep- 
arable from the others iff AjcrAj + ^^ > 0, where Aj 
is a diagonal matrix implementing the partial transpo- 
sition at the level of CM, i.e. inverting the sign of the 
j-th momentum. In our case Ai = Diag(l, —1, 1, 1, 1, 1), 
A2 = Diag(l, 1,1, -1,1,1), A3 - Diag(l, 1,1, 1,1,-1). 
When we have a non Gaussian state the violation of the 
above condition is still a sufficient condition for entan- 
glement, whereas the necessary part is lost. 

Let us start by focusing attention on a A'^-atom pure 
subradiant state ^. The density operator is gp — 
\sr)pp{sr\ and the only non- vanishing terms of the CM 
(T[gp] are the diagonal ones, which are reported in ap- 
pendix |Al This fact is a direct consequence of two se- 
lection rules affecting the geometrical form of subradiant 
states: the first one is the conservation of the number 
of atoms and the second one is the presence of the k- 
index in every Fock component of \sr)p. As a matter of 
fact, the CM (T[gp] satisfies the three separability con- 



ditions AjCrAj + |J7 > 0, j = 1,2,3 and therefore we 
have no information about the entanglement properties 
of pp, i.e. we cannot use the PPT condition translated to 
continuous-variable systems to analyze and quantify en- 
tanglement. On the other hand, we know from the pre- 
vious analysis that gp is genuinely tripartite entangled 
state and this suggests that it should be far from being 
a Gaussian state. Indeed, the characteristic function is 
given by a Gaussian modulated by Laguerre polynomials 
and thus shows a distinctive nonGaussian shape. 

In the following we will evaluate quantitatively the 
nonGaussianity of the subradiant states upon the use of 
the nonGaussianity (nonC) measure introduced in [25|, 
which is defined as 



m 



Dhs[Q^A _ /^e + Mr-2Kg 



Me 



2Me 



(15) 



where DhsIQtt] denote Hilbert-Schmidt distance be- 
tween the state g under scrutiny and its reference Gaus- 
sian state T, i.e. a Gaussian state with the same mean 
value X[t] = X[g\ and CM (t[t] = (T[g\, fig = Trlg'^] is 
the purity of g, and k^t — Tr[£iT] is the overlap between 
g and r. The nonG measure 6[g\ possesses all the prop- 
erties for a good measure of the non-Gaussian character 
and, in particular, 6[g] — iS g is a. Gaussian state. For 
the pure subradiant states gp = \sr)pp{sr\ the reference 
Gaussian state has the form t = vq iSi vi 'Si V2 where Vj = 
(1 - yj)Y.syj\^)i'^\^ Vj = ^i/(l + ^j) is a (Gaussian) 
thermal state with Nj average thermal photons. For a 
three-mode thermal state we have (111 = (722 = |(A'o + ^), 
CT33 = tT44 = i(A^i-l-i), 0-55 = 0-66 = ^iN2 + ^). Since for 
a subradiant state the average occupation numbers are 
given by A^o = p- {k)p, Ni = 2{k)p, N2 = N-p- {k)p we 
have, after some calculations, that the non-Gaussianity 
of gp is given by 



2 2 



Vj 



k=0 



(16) 



For a large number of atoms N '^ 1 the nonGaussianity 
of ^[gip] may be evaluated upon exploiting the Eq. ^ pH)) 
i.e. the relation between p and e. In practice, every value 
of e corresponds to a single p and in turn to the pure 
subradiant state gp. In Fig]?] (top) we show the behavior 
of the nonGaussianity S[gp] for A^ = 50 and < e < 1 -I- 
v2. As it is apparent from the plot, the subradiant states 
are always non-Gaussian and the value of S is almost 
constant. For A^ = 2, 3 it is possible to evaluate explicitly 
the non-gaussianity for the stationary state (JH). Since 
the number of atoms is small p and e are independent 
one another and we are not allowed to use relations (jH])- 
(ITUl) . For A^ = 2,3 we have the ground state, go, and a 
single subradiant state, gi, and the nonGaussianity may 



be expressed as 
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The reference Gaussian state is again a three-mode ther- 
mal state, as it can be easily demonstrated by noticing 
that the CM of Qs is the convex combination of those 
corresponding to ga and gi i.e. (T[gs] = Picr[pi] + (1 — 
Pi)cr[£io] = Pi- Pi is given by eq.Q for A^ = 2 and it 
may be evaluated numerically for A'^ = 3. The overlaps 
are given in appendix |B] The behaviour of the nonCaus- 
sianity d[gs] as a function of e for N = 2,3 is shown 
in Fig. m (bottom). Notice that the behaviour of S[gs] 
does not show a strict correlation to that of the negative 
eigenvalue A{e) and thus we cannot use nonCaussianity 
to assess quantitatively the failure of CV PPT condition 
in detecting stationary entanglement in our system. 



IV. CONCLUSIONS 



Appendix A: CM of the subradiant state 

The covariance matrix of the subradiant state gp = 
\sr)pp{sr\ is diagonal. The nonzero matrix elements may 
be expressed as 



C^ll = (722 = 



k=a ^ ' 



(Al) 



\{v-^k)v-\ 



a33 = a44 = ^E/3^(4fc + l) 



fe=0 



(fc)p + 4 



1 



-^E/3^f^-p^fc + ^ 



k=a 



(A2) 



(A3) 



iV-p-(A;)p + - 



where we have defined 



In this paper we have considered a systems made of A^ 
two-level ultracold bosonic atoms in a ring cavity where 
subradiance appears in the bad cavity limit. We pro- 
vided a full quantum description of the dynamics and 
have shown the appearance of stationary entanglement 
among atoms. We have evaluated the amount of steady 
state entanglement for N = 2 and A^ = 3 atoms as quan- 
tified by the negativity of the partially transposed den- 
sity matrix. We have also investigated the entanglement 
properties of the systems as a continuous variable one 
and have shown that it is not possible to detect entan- 
glement due to the nonCaussian character of subradiant 
states. It is important to notice that the life time of 
experimentally demonstrated entangled states is gener- 
ally limited, due to their fragility under decoherence and 
dissipation. Therefore, in order to decrease dissipation 
and decoherence, strict isolation from the environment is 
usually considered. On the contrary, in our system, the 
action of the coupling to the environment drives it into 
a steady entangled state, which is robust with respect to 
the interaction parameters. Our results pave the way for 
entanglement characterization [2^ of subradiant states 
and suggest further investigations for large number of 
atoms. 




{2k)\{N ~ p - k)\ 



{p~k)\ 



and (fc)p = |Cpp X)fc=o Pk^ ^^ ^^^ mean value of k for the 
p-subradiant state. 



Appendix B: Overlaps 

The overlaps Kg-r between the states gj and their ref- 
erence Gaussian states for A" atoms, which appears in 
the expression of the nonGaussianity 5[gs\ are given by 



^eo,T - '^Qo.r - 2/2 



Hi 



1-^ 



c(2) = 



c(3) = 



yl + 2e^yo 2/2 -pr 1 - Vj 
1 + 2e2 11 1 + y^ 



i=o 



y2 {yl + 4e2yo yz) T-r 1 - % 



l-h4e2 



nf 



3=0 



■Vj 
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